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Abstract 


User equilibrium is a central concept for studying transportation networks, and one 
can view it as the result of a dynamical process of drivers’ route choice behavior. In 
this paper, based on a definition of 0-D First-In-First-Out violation, we propose a new 
dynamical system model of the route choice behavior at the aggregate, route level for 
both static and dynamic transportation networks. An equilibrium of such a dynamical 
system can be a user equilibrium or a partial user equilibrium. We prove that, for 
static, symmetric traffic assignment problem with fixed or variable demand, only user 


equilibria are stable for the dynamic al system, and the objective function in t’ 


ematical programming formulation (jBeckmann. McGuire, and Winstenl . 


1956 


le math- 
) can be 


considered as the potential energy of the dynamical system. We then present an Euler- 
based perturbation method for finding user equilibrium and solve two examples for 
both static and dynamic traffic assignment problems. This new model is simple in 
form and could be applied to analyze other properties of transportation networks. 
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1 Introduction 

The traffic assignment problem (TAP) is to assign origin-destination (0-D) flows to a road 
network. Solving this problem is a core step of four-step transportation planning method 


and essential for analyzing congestion and other performance of a networ k. It is general 


believed that assigned traffic flows are in the state of user equilibrium (UE) fjWardropl. Il952l) . 


y 


where “the journey times of all routes actually used are equal and less than those which would 
be experienced by a single vehicle on any unused route”. That is, for the same 0-D pair, 
all used routes share the same travel time, and unused ones are not shorter. A fundamental 
task in solving TAP is to compute travel times on all links. In static TAP, it is assumed that 
there exist so-called link performance functions, which map link flows to travel times. For 
static TAP, link performance functions can be separable, symmetric, or monotone. While 
in dynamic TAP, not only the number of vehicles but the dynamics of traffic can affect the 
travel time on a link, and there may not exist any meaningful link performance functions 


flDaganzol . 


1 995bh 


There are four well-known formulations of TAP: namely, mathematical pr ogramming 


varia tional inequality, nonlinear complementarity, and fixed point formulations (jPatrikssonl . 
I 994 J ). In addition, there exists another so-called dynamical system approach, inclu ding dy¬ 


namical models of stochastic, individual route choice behavior ( gg. 


namical mo dels of aggregate, deterministic change of route flows (jSmithl . 


Horowita. 


1996 


1984; 


1984) and dy- 


Nagurnev and Zhand . 


19971 ). Equilibria of these dynamical systems are usually UE. However, these dynami¬ 


cal systems are very complex in form and have not been widely used in solving and analyzing 
static or dynamic UE for large-scale transportation networks. 
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In this paper, we propose a new dynamical system formulation of TAP. For static, sym¬ 
metric link performance functions, we define a new aggregate, deterministic route choice 
dynamics based on a so-called First-In-First-Out (FIFO) violation function among different 
routes of the same 0-D pair. We then show that, in equilibria of this dynamical system, 
“the journey times of all routes actually used are equal”, but may not be “less than those 
which would be experienced by a single vehicle on any unused route”. That is, an equilib¬ 
rium of this dynamical system may be a UE or a partial user equilibrium (PUE). For static, 
symmetric TAP, however, we are able to show that PUE are always unstable, and only UE 
are stable. Then, based on the stability properties of UE and PUE, we propose a numerical 
algorithm for computing stable UE. We further extend this new dynamical system formula¬ 
tion for TAP with variable demand and for dynamic TAP that incorporates traffic dynamics. 
Through this study, we intend to develop a model that can be used for theoretical analysis 
of properties of various transportation networks. 

The rest of the paper is organized as follows. In Section 2, we introduce a new dynamical 
system of static TAP and analyze stability properties of its equilibria. In Section 3, we present 
a hnite difference method and a perturbation-based method for computing equilibria of the 
dynamical system and UE. In Section 4, we extend the new dynamical system formulation 
for dynamic TAP. We make some concluding remarks in Section 5. 


2 A dynamical system model of static traffic assign¬ 
ment problem 


For static 
to that in 


AP, we ado pt the network notation system given in Table [H which is similar 


(jSheffij, 


1984J ). but with several new dehnitions. We also have the following basic 
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relationships: 


and 


Qrs = ^ r,S 

h 

(la) 

K 

/[" >0, V fc, r, s; 

(lb) 

= E E/rc*. 

(2a) 

r,s k 


= Vi.r.s, 

(2b) 


We call f as a UE state in the sense of flWardropl. 


1952 ) if and only if for f G JF 


:) = 0, V fc, r, s 

(3a) 

> Vrs, V fc, r, s. 

(3b) 


As w e know, UE is the solntion of the following mathematic programming problem fiReckmann. McGuire, and 
19561) 


mm 2 ; X = 


/ 

/ ta{uj)du, 


(4) 


whose objective fnnction is hereafter called the BMW objective fnnction. The BMW ob- 


i ective fnnction 2 :(x) is strict 


(jPafermos and Sparrowl . 


y convex in x for static, symmetric TAP with hxed demand 


19691 ) ■ In this section, we assnme that z(x) is convex, bnt not nec¬ 


essarily strictly convex. That is, z(x) still attains its minimnm at UE solntions of link flows, 
bnt UE solntions may not be nniqne in x. The BMW objective fnnction can also be written 
in ronte flows, z{f) = 2 ;(x(f)), and, at UE, z{f) also attains its minimnm. 


2.1 O-D FIFO violation and a dynamical system model of aggre¬ 
gate route choice behavior 

On a road link, EIEO among vehicles is violated when one vehicle passes another. Snch link 
FIEO violation can be cansed by heterogeneity between vehicles, different traffic conditions 
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on different lanes, or traffic signals. M easurements were deve l oped for link FIFO violation 


among diffe rent groups of vehicles in fjjin and Javakrishnanl . 


vehicles in fj.Tin et al. 


200511 and among individual 


20061). For vehicles of the same 0-D, even when vehicles observe the 


FIFO principle on all links and routes, the FIFO principle can still be violated if vehicles 
using different routes experience different travel times. Here we are interested in such 0-D 
FIFO violation and assume that in 0-D FIFO state vehicles departing from the origin at 
the same time arrive at the destination at the same time. Hereafter, FIFO violation means 
0-D FIFO violation if not otherwise noted. 

To show how to measure 0-D FIFO violation in static transportation networks, we con¬ 
sider a network with one 0-D pair, r — s, and two alternative routes. In one assignment, 
one obtains two differen t route travel times, The cumula tive arrival 


curves at the origin and destination (iMoskowitz and Newmanl . 


1963; 


Newell, 


19931) are shown 


in Figure [H where Vrs is the average 0-D travel time. Then we can define FIFO violation 
for two route flows, by the areas of the two regions shown in Figured! but subject 

to opposite signs. Note that here the unit of FIFO violation is vehicle x time. 

For a general road network, we can define the FIFO violation function for route k con¬ 
necting 0-D pair r — s by (V fc, r, s) 


r* — n 


rsfricl^-Vrs), 


(5a) 


I 


whose unit is vehicle^ x time □. Since the average 0-D travel time is 


Vrs — 


E frs^rs 

j ^ j ^3 


(5b) 


the FIFO violation function can be re-written as 


jrs \ ^ frs rs\ rrs \ ^ nrsf rs 

Jk / j j j ) Jk / J J j \^k )’ 


( 6 ) 


^Here we multiply FIFO violation defined in Figured] by for convenience, and our discussions apply 
if all 0-D FIFO violation functions of an 0-D pair are multiplied by the same constant. 
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To model the route choice dynamics at the aggregate, route level, we introduce a new 
principle so that the traffic flow on a route is decreased when the travel time on the route is 
longer than the average 0-D travel time. In particular, assuming that all route flows depend 
on a decision variable r, whose unit is {vehicle x time)~^, we propose that the rate of change 
in a route flow be negative to the corresponding 0-D FIFO violation; that is, 

- /r = J" = fk - y)/;- = q„K‘(ci- - t,„), -i k,T,s (7) 

j 

where is the derivative of with respect to r. From (JTj), we can have the following 
autonomous dynamical system of route flows f 

-f = J(f). (8) 

Note that, with the definition of 0-D FIFO violation in flSal) . we have a negative sign on the 
left-hand side of (IH]). 

The new dynamical system, (0), has the following properties, (i) If fl^{0) = 0, then 
— 0 "T > 0; that is, if a route is unused initially, it is unused at any decision 

step, (ii) The dynamical system is feasible in the following senses. First, if initial route flows 
are non-negative, they are always non-negative at any r, since the trajectories of 
cannot cross the boundary /" = 0. Second, 0-D flows are conserved; i.e., = ^rsi 

since Xlfc/"(''’) ^ ^ 

Definition 2.1 (FIFO Equilibrinm) Roots of equations J(f) = 0 are equilibria of the 

dynamical system Since 0-D FIFO violation functions equal zero in these equilibria, an 
equilibrium of the dynamical system is also called a FIFO equilibrium. 

Remarks. We have the following properties of a FIFO equilibrium. 

1. f G JF is an equilibrium of ([8]) if and only if it satisfies one of the following conditions: 
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(a) Jr = rr EM' - C^)/^ = O (Wk, r, s), 

(b) Jr = Qrs^r ~ 'J/)^r ^ ^ i'Jk,r,s), where ^r = fr/Qrs is the proportion 


of flow on route k among flow on 0-D pair r — s. 



,rs 


k 


Vrs) = 0 (Vfc,r, s), 




< •••, where fr > 0,V fc 


2. The dynamic system usually has multiple equilibria. For example, any vertex in the 
convex polygon dehned by ([T]) is an equilibrium. In equilibria, used routes share the 
same travel time, but unused routes may have longer or shorter travel times. 

3. It is straightforward that, if f is a UE, then it is a FIFO equilibrium. However, some 
equilibria may not be UE. We call such FIFO equilibrium as partial user equilibrium 
(PUE), since it can be considered as UE for a partial set of routes. We can see that a 
PUE satishes the hrst part of Wardrop’s dehnition of UE; i.e., in a PUE, “the journey 
times of all routes actually used are equal”. In real-world transportation networks, 
such PUE states could be realized if drivers prefer familiar routes over shorter ones 
connecting the same 0-D pair. 


2.2 Stability properties of UE and PUE 


In the following, we show that UE and PUE have different stability properties. First, we 
prove the instability of PUE. 

Theorem 2.2 A PUE is unstable for the dynamical system 
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Proof. If f is a PUE, there exists an unused, shorter route k for some 0-D pair r — s. Then, 
for a used route j, c]f < = 0, and > 0. We introduce a perturbation to this 

equilibrium by shifting cq (/J^ > cq > 0) from route j to route k. We denote the resultant 
state by f and the corresponding vector of route travel times by c. Since route travel times 
are continuous in f, we can have a sufficiently small Cq such that route k is still shorter than 
other used routes; i.e., < dy and d'ff < c[® for any f[^ > 0. Then at = eo, the local 

changing direction in is 

/r = >0- 

i 

That is, the perturbed state f will drift away from the original equilibrium, and the original 
PUE is not stable. ■ 

Then, we introduce a Lyapunov function for the dynamical system ([8]) as follows. 

Theorem 2.3 The function 


V (f) = z{i) — min z{i) 


(9) 


is a Lyapunov function ALaSalle 
is 


JML; 


Stroaatz 


1994 ^ of the dynamical system at UE. That 


1. U(f) > 0, V f G JP and f is not a UE; 

2. U(f) = 0 if and only i/f G P, where £ is the set of UE; 

3. -grad U(f) ■ J(f) < 0 z/f is not a UE. 

This is equivalent to saying that the BMW objective function can be considered the potential 
energy of the dynamical system 

Proof. From the dehnitions of the Lyapunov function V (f) in ([2]) and the BMW objective 
function z(f) in (jl]), we can see that the hrst two statements are correct. We prove the third 
statement as follows. 
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We first compute the gradient of V (f) with respect to f, whose (rs, /c)th element is 
(V k,r,s) 


Therefore, 


^V{^) 


dz{yi) 

'aJjF 




—grad V (f) ■ J(f) 


which leads to 


\ ^ \ ^ ^ ^ rrs _ \ ^ \ ^ rs rrs \ ^ / rs 

/ V / V Qfrs ~ ~ / , / , Ik / , \^k )Jj ) 

rs k ^ ^ rs k j 


-gradV{f)-m = - EE EW “< 0 , ( 10 ) 

rs k j>k 

since all route flows are non-negative. From the dehnition of equilibria, we can see that 
—grad ld(f) • J(f) = 0 if and only if f is an equilibrium, and —grad ld(f) ■ J(f) < 0 for 
non-equilibrium states. Therefore, V (f) is a Lyapunov function of the dynamical system ([H]). 


We have the following property for the set of equilibria of the dynamical system (j8|). 

Lemma 1 The set of equilibria is closed. That is, given a sequence fj —>• P, if all fj are 
equilibria, then f* is also an equilibrium. 

Proof. Since (V k, r, s) 

4'(f) = Eo«(EE'*aT/r) 

a mn I 

is continuous in f, FIFO violation functions, J(f), are also continuous in f. Then, with 
fj —f*, J(fj) = 0 implies J(f*) = 0. That is, f* is also an equilibrium, and the set of 
equilibria is closed. However, note that the set of equilibria is generally not connected, since 
we may have multiple isolated regions of equilibria. ■ 

Further, we have the following properties for the set of UE. 

Lemma 2 The set of UE, S, is closed and connected. 


9 





Proof. Given a sequence of f* G and fj —> f*, then from continuity of the BMW objective 

function, we can see that z(fj) —> ^(f*)- Therefore, is also minimum, and f* G S. Thus 
S is closed. Since z{x.) is a convex function and x is in a convex set, the set of UE, S, is 
connected in x. Further, since Xa = Xlrs '^k ^a\fk^ ^ continuous mapping from T to A, 
we have that S is also connected in f. ■ 

Finally, we show that UE is stable for symmetric link performance functions as follows. 

Theorem 2.4 A UE is locally, asymptotically stable for the dynamical system ^ in the 
sense that solutions of ^ converge to a UE for a non-UE initial state that is close to the 
UE. 


Proof. From the lemmas above, we can End a region around UE, where all non-UE states 


are not equilibria. Then —grad U(f) ■ J(f) <0 for all the non-eq uilibriuin sta 


;es in this re¬ 


gion. Therefore, according to the Lyapunov’s stability theorem in (iSmithl. I1984J ). asymptotic 
solutions of (IHl) converge to UE solutions for non-UE initial states. That is, the set of UE 
is asymptotically stable. Note that we can have multiple UE and cannot use the traditional 
Lyapunov’s stability test. ■ 

From Theorems 12.21 and 12.41 we can see that UE are the only stable equilibria of the new 
dynamical system for symmetric link performance functions. Note th at, however, we can 


have unstable UE when link performance functions are non-monotone (INetter . 


1972h . 


2.3 A dynamical system model of static TAP with variable de¬ 
mand 


The dynamical system (0) for fixed demand can be extended for st atic TAP with variable 


demand, whose mathematical programming formulation is given by flBeckmann et ah 


1956 
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Sheffil . 


1984h 


^Qrs 


min^ 


Xq) = 


(11a) 


subject to ([I]) and 


> 0, V r, s, (11b) 

where Urs{qrs) is a decreasing travel time function associated with 0-D pair r — s and 

the inverse of the demand function. For symmetric link performance functions, the BMW 
objective function in filial) is still convex, and there exists a unique user equilibrium in x 
and q. 

For this problem, we introduce the following dynamical system 

/r = -qrsfricl'-Urs), (12a) 

qrs = -qrsC^fr^l'-qrsUrs), 'ir,S. (12b) 

k 

Or, if introducing the average 0-D travel time v^-s as in fl5bj) . we can have an equivalent 
dynamical system 

/r = (13a) 

j 

qrs = -qrsi'l^rs - Urs), '^r,S. (13b) 

Similar to ([8]), this dynamical system can have both UE and PUE as its equilibria. Further, 
we can have that: (i) PUE are unstable; (ii) the BMW objective function in (IIlap can be 
considered as the potential energy function of the dynamical system ca; and (hi) only 
user equilibria are stable. These statements can be proved in the similar fashions as in the 
preceding subsection. 


11 









3 Computation of static UE 

In the preceding section, we have derived a physically meaningfnl model (IH]) for static TAP 
and analyzed its stability properties. In this section, we introdnce a nnmerical method for 
compnting UE based on the dynamical system formnlation. The pnrpose here is to show 
that the new formnlation is feasible for finding UE nnmerically. Therefore, we do not concern 
with the compntational efficiency of this method, which is an important factor for solving 
large-scale problems bnt not essential for analyzing median-sized prototype networks. 


3.1 Algorithms 


We first discuss algorithms for solving UE of static TAP with fixed demand. Since an unused 
route in the initial state by vehicles of an 0-D pair is always unused, we only consider those 
paths with non-zero initial flows and denote the number of initially used routes connecting 
0-D pair r — s by Krs (V r, s). 

By discretizing the dynamical system ([8]) in the decision space as 


f(r -I- At) — f(r) 
At 




we can update f by 


f(r -h Ar) = f(r) — ArJ(f(r)), 


(14) 


whic 


1 is the standard Euler’s method for solving ordinary differential equations fjStrogata. 


199411 


Note that, at each step, we need to re-assign intermediate route flows to links to compute 
route travel times and obtain FIFO violation functions J(f). Therefore, the total computa¬ 
tional load of the finite difference method is proportional to the number of steps and that 
of computing FIFO violation functions. The index of convergence can be simply defined as 
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the norm of FIFO violation functions as follows 




( 15 ) 


which equals to zero if and only if f is an equilibrium. 

For an arbitrary initial state f(0), the equilibrium of the dynamical system ([H]) may not 
be UE, but a PUE, in which we can End unused shorter routes. For a PUE, we can perturb 
it by shifting proportions of flow from a used route to the shorter unused routes and obtain 
another equilibrium, since a PUE is unstable in these perturbation directions, as shown in 
the proof of Theorem 12.21 We repeat this process until we reach a UE, where there are 
no shorter unused routes. The flow-chart of the algorithm is given in Table [2l Note that 
the initi al guess does not have to be an all-or-nothing assignment as in Frank-Wolfe’s (FW) 


method fISheffil . 


1984) . Rather, with iP-shortest routes as initial state, we could find UE with 


few or no perturbations. 

For static TAP with variable demand, corresponding to the two formulations (|T^ and 
fin]) we can have two types of solution methods. In one, with an initial guess of qrs (Vr, s) and 
fr we directly solve fll2al) and fll2bD . by using the aforementioned perturbation- 

based method. In the other, with an initial guess of qrs, we solve the corresponding UE 
for a fixed demand problem (11 dab with Vrs = Urs, and then update qrs and with the 
difference between Vrs and Urs by fll3ap and fll3bp respectively. 

3.2 An example 

In thi s subs ection, we study TAP on a simple network given by (Fig. 5.1 on page 114 of 


Sheffil . 


1984) . which is shown in Figure [2l For this network, FIFO violation functions can be 


written as 


Ji = Xi ((lO(l + 0.15(y)") - 20(1 + 0.15(^)")) 0(2 + (lO(l + 0.15(y)") - 25(1 + 0.15(y)' 
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J 2 = X 2 (( 20(1 + 0.15(^)") - 10(1 + 0.15(^)^)) + ( 20(1 + 0.15(^)") - 25(1 + 0.15(^)^)) Xg) 

J 3 = X 3 (( 25(1 + 0.15(^)") - 10(1 + 0.15(^)")) 0:1 + ( 25(1 + 0.15(^)") - 20(1 + 0.15(^)")) x^) 

and the dynamical system of the network is (i = 1, 2, 3) 

Xj t/j. 

First, we consider an initial state at (xi, X 2 , Xg) = (3.39, 5.00,1.61) and the corresponding 
route costs (ci, C 2 , cg) = (22.3, 27.3, 35.3). In the region of X 1 +X 2 < 10 (xg = 10—xi—X 2 ), the 
convergence is shown in Figure[3], and the solution trajectories in Figure IHwith r = 0.02 and 
Ar =0.0005, 0.001, and 0.002, respectively. From these hgures, we can see that convergence 
rates and solution trajectories are similar for different decision step sizes, as long as the hnite 
difference form of ([8]) is stable (i.e. Ar is sufficiently small). In addition, convergence rates 
are constant when solutions get closer to equilibria, and there is no zigzagging effect. 

With different initial conditions, we can find all seven equilibria for this network as 
shown in Table |3l We perturb the first six PUE by shifting 0.05 flow from a used route 
to each unused shorter route. After perturbation, we use the finite difference method with 
Ar = 0.0005 and r = 0.1 to obtain solution trajectories of the dynamical system ([H]). 

As shown in Figure [5l all solution trajectories converge to UE after perturbation. This 
example confirms the theory that only UE are stable and demonstrates the feasibility of the 
perturbation-based method for finding UE. 

4 A dynamical system model of dynamic traffic assign¬ 
ment problem 

In this section, we consider dynamic TAP, in which all quantities are time-dependent. Since 
the concepts and definitions are all similar to those in static case, we do not specifically 
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include “dynamic” in them and keep the same terminologies as before. For example, UE 
is meant to be dynamic UE. We adopt the network notation system given in Table 0] and 
have the following remarks. First, flow-rates are derivatives of flows, and their units are 
number of vehicle per unit time and number of vehicles respectively. Second, all variables 
depend on the decision variable r, although not explicitly indicated. Third, the formulation 
is route-based, and we do not use link flows, link travel times, or indicator variables. In 
addition, we have the following basic relationships: 


Prs{r, t) = Y^ /r(u t), Prsis, t) = Yl ^), V r, s 


/r(u t) > 0, /["(s, f) > 0, V k, r, s; 


V r,s 

(16a) 

V k,r,s-, 

(16b) 


where Vm and ar e cumulative arrival curves at origins and destinations (iMoskowitz and Newmanl . 


19631: 


Newell 


1993 ). and 


grs{r,t) = 

(r,t), qrs{s,t) = ^gl^{s,t), V r,s 

(16c) 

K 

K 

gl"{r,t) > 0, gl"{s,t) >0, V k,r,s. 

(16d) 


4.1 Formulation 

For dynamic TAP, we propose the following dynamical system for adjusting route flow-rates 

V k,r,s (17) 

where the time-dependent 0-D FIFO violation function is dehned by 

= grs(r,t)gj:''(r,t)(4"(r,t) -Vrs(r,t)), Vk,r,s (18a) 

and the 0-D average travel time for vehicles departing origin r at time t, Vrs(r,t), is 

Vrs(r,t) = ^gj"(r,t)cj"(r,t)/qrs(r,t). (18b) 

j 
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Note that fflTj) is equivalent to — V = — t). This dynamical system, continuous in 


time, is a direct extension of ([8]) at any time t. Althongh there nnlikely exist link 


fnnctions when considering capacity constraints and link interactions (jPaganzol. 


perform ance 


1995hh . the 


ronte travel time c^®(r, f) can be compnted from a model of network vehicnlar traffic. For 


exa mp 


in fjjin 


e, we can use the commodity- based kinematic wave model of network traffic developed 


2003; 


Jin and Zhand. 


200411 ■ in which FI 


^0 violation among vt 

Tides 

Jin and Javakrishnan. 

2005) 


If 0-D FIFO violation functions equal zero for all r, s, k, t, the corresponding state of flow- 
rates gl.^{r,t) (V k,r,s) is called a dynamic equilibrium. We can have multiple equilibria, 
and trivial equilibria include those in which only one route is used for an 0-D pair. Those 
equilibria are called dynamic partial user equilibria (PUE) if there are unused shorter routes, 
and dynamic user equilibria (UE) otherwise. This dehnition of dynamic UE is consistent 
with Wardrop’s. Further, from properties of static user equilibria, we have the following 
conjecture: Dynamic PUE are unstable for (1T7|1 : i.e., stable equilibria of (|T7|1 are dynamic 
UE. In general, we expect to have multiple dynamic UE, and a dynamic UE can be unstable, 
due to non-monotonicity of link travel times with complicated interactions between traffic 
streams in a dynamic road network. 


4.2 Computation of dynamic UE 

In this subsection, we describe a computational method for Ending equilibria and stable 
UE based on the new formulation dzi. Here we assume that origin demands are given 
during an assignment time duration [0,To] and zero outside the interval. We first discretize 
the assignment duration into n time intervals with time instants = nAt, where n = 
0,1, 2, ■ • •, iV and N = Tq/A t. We assume that during [f„, the in-fiow-rate of route k 
connecting 0-D pair r — s is constant, gl.^{r, tn), and in-flows, fl^{r, tn), are piece-wise linear 
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functions. 

In order to compnte ronte travel times, we can nse any traffic simnlator with traffic 
simnlation time dnration of [0, T] (T > Tq), dnring which all vehicles shonld be able to hnish 
their trips. Also, in order to control nnmerical errors when compnting travel times, we nse 
smaller simnlation time steps, At/M (M > 1), and simnlation time instants of 
{i = 0,1,2, ■■■, MNT/Tq). Then, given an initial gness of g/f{r,tn) (V k,r,s), bonndary 
conditions at destinations, and initial traffic conditions, we can nse a traffic simnlator to 
obtain solntions of arrival cnrves at destinations f/,^{s,t) at any simnlation time instant. To 
compnte c]f{r,tn)-i which is the average travel time of vehicles from fjl^{r,tn) to fl^{r,tn+i), 
we hrst hnd time instants tm„ and snch that tm^) and are the closest 

to fl^{r,tn) and fl^{r,tn+i) respectively. As illnstrated in Fignre[6l the total travel time of 
vehicles departing dnring {nAt, {n + l)At) is the area bonnded by the two cnmnlative flow 
cnrves at the origin and destination r — s and can be computed by 


cl"ir, tn)gl"{r, tn)At = 


r/r(utn+i) 

hr(rM) 


« (/r(r,(„+i) - -n- 1)A« 




i=mn + l 


from which we can compnte the approximate 0-D FIFO violation J/f{r,tn) by nsing flTSll . 
Then we solve fll7p by a hnite difference approximation as follows 


£/r(A^n)|r+Ar = gk{r,tn)l 


- 4" (a ^n)Ar, 


V k,r,s. 


( 20 ) 


where n = 0,1, ■ • A — 1. From these solntions, we can constrnct a new time series of 
in-flow-rates, gl/{r,tn), a.t t + At. Here we nse the following the index to measnre the 
convergence of solntions 


Er.EfcEn=l('/r(A tn)) 


NE 


rs ^rs 


(21) 
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where Krs is the number of initially used routes connecting 0-D pair r — s. If the hnite 
difference equations are stable and convergent for a sufficiently small At, solutions will 
converge to equilibria for a sufficiently large r. The flow-chart of this algorithm for finding 
an equilibrium is given in Table [5l Further, if there are shorter unused routes in a PUE, we 
shift some flows to such routes and find another equilibrium, until we reach a stable UE. 


4.3 An example 


We study dynamic TAP in the n etwork shown i n Figure [71 where t wo routes have the same 


triangular fundamental diagram fiMnnial et ah 


1971 


Newell, 


19931) . The network is empty 


initially, and there is no capacity constraint at the destination. Here we consider a constant 
arrival flow-rate at the origin, q^, during time interval [0,To]. That is, qrs{r,t) = go for 
t G [0, To] and 0 otherwise. 


We use the Lighthill-Whitham-Richards (ILighthill and WhithanJ . 


1955; 


Richards! . 


195d) 


traffic model to analyze traffic dynamics on each link. For example, for link 1, however many 
vehicles are waiting to enter, the maximum flow-rate is gc = 1. Since there is no bottleneck 
on the link or the destination, we know that traffic on link 1 is free flow, and the travel 
time is 1. Although the link travel time is constant, the waiting time at the origin is time- 
dependent. If the arrival flow at the destination for route 1 is fi{s,t), then the departure 
flow at the origin is /i(s, 1 +t). Denoting = iAt/M, we then have the following equation 

fi{s,l + ti+i) = fi{s, l + U) + ^mm{qc, + + gi{r,ti)},{22a) 

where /i(s, 1) = 0, /i(r, 0) = 0, and _j_ jg the maximum flow-rate 

that can be sent from the origin if there is no capacity constraint on this link. This method o 


computing boundary fluxes is based on t 


1995a : 


Lebacqnel . 


19961 : 


Jin and Zhangl . 


le su pply-demand method developed in flDaganzo 


20031 ). Thus, (I22ap gives the dynamic model for 


finding arrival flows at the destination on route 1. Similarly, we can have the following 
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model for route 2: 


/2(s,2 + ti+i) = f2{s,2 + ti) + ^max{q^ 


f2{,r,U) - f2{s, 2 +U) 
At 


M + g2{r, ti)}l22h) 


where / 2 (s, 2) = 0 and / 2 (r, 0) = 0. 

We set To = 1, T = 8, iV = 20, M = 10, qo = 5, t = 160, and At = 0.05, and consider 
deterministic initial time series of gi{r,tn) = ego and g2{r,tn) = (1 — c)go (n = l,---,iV). 
When c = 0.5; i.e., both routes have the same flows initially, solutions of cumulative arrival 
curves of route flows and total flows at the origin and the destination are shown in Figure 
[HI from which we can clearly see that these solutions satisfy the 0-D FIFO principle. At 
different time instants, the travel times are shown in Figure [9], from which conhrms that the 
solution is an equilibrium given by 




5t, 0<t< 0.25; 

< 1.25+ 2.5(t- 0.25), 0.25 <t<l; 
3.125, t>l, 


(23a) 




f2{r, t) 


0, 0 < t < 0.25; 

< 2.5(t- 0.25), 0.25 <t<l; 
1.875, t > 1. 

V ’ 


(23b) 


Since there are no unused shorter route in the solution, the equilibrium in (123|1 is a UE. 

With the aforementioned deterministic initial time series with c = 0.5, 0.95, and 0.05, and 
another random initial time series, we hnd that all solutions converge to fl2^ and the index 
of convergence is shown in Figure fTOl from which we can see that solutions converge slowly 
at the beginning stage, then exponentially when it is close to fl23l) . Note that gi(r, t„) = go 
(V n = 1, ■ ■ ■, N) and g 2 {r, tn) = qo (V n = 1, ■ ■ •, A^) are two trivial equilibria. Since initial 
time series with c = 0.95 can be considered as a perturbation around gi(r, t) = go and those 
with c = 0.05 a perturbation around g 2 {r,t) = qo, we can conclude that the UE in fl23|) is 
stable. Convergence with random initial time series further confirms the stability of the UE. 
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5 Discussions 


In this study, based on a definition of 0-D FIFO violation among different routes connecting 
the same 0-D pair, we derived and analyzed a new dynamical system formulation of the 
traffic assignment problem for both static and dynamic road networks. For static, symmetric 
link performance functions, we also showed that the dynamical system is stable at UE, 
and the well-known BMW objective function can be considered the potential energy of the 
new dynamical system. Through this study, we can see that the new model is physically 
meaningful, theoretically rigorous, and numerically feasible. 

The new dynamical system model is simpler in form than existing aggrega t e, det ermin- 


istic dynamical systems proposed in fjSmithl . 


1984 : 


Nagurnev and Zhand . 


1996, 


3 ). 


199711. Also 


different from existing dynamical systems, its equilibria can be UE or PUE. Compared with 
mathematical programming formulations, this formulation is applicable for non-monotone 
traffic assignment problem, and we can apply this dynamical system model to analyze the 
stability of multiple UE. In addition, this new approach can be extended for other traf¬ 
fic assignment problems. For example, since system-optimal assignments are equivalent to 
user-equilibrium assignments with modified link performance functions, we can also have a 
similar dynamical system formulation of system-optimal assignments. In addition to being 
applied to theoretical formulation and analysis of various traffic assignment problems, the 
new model can also be used to numerically study small- to medium-sized networks, which 
can be solved in a reasonable amount of time with modern computers. 

In the future, we will be interested in further theoretical and numerical investigations of 
the new dynamical system. Theoretically, the existence of stable equilibria of the dynamical 
system of static TAP with monotone link performance functions is subject to further inves¬ 
tigations; it is also possible to incorporate network control policies and drivers’ familiarity 
and preference of certain routes into the dynamical system; it will also be interesting to 
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study the existence, uniqueness, and stability properties of dynamic UE in a general net¬ 
work with this formulation. Numerically, although the Euler-based perturbation method 
proposed in this study is shown to be feasible with two simple examples, the method only 
converges linearly and is not sufficiently efficient for large-scale problems. In the future, we 
will be interested in developing more efficient methods based on this formulation. For exam¬ 
ples, we could develop parallel computational methods, or hybrid methods integrating the 
Frank-Wolfe method and the perturbation-based method; we can also use the fourth-order 
Runge-Kutta method, which is more efficient than the Euler’s method. 
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M 

A 

7^ 

5 

/Crs 

T 

At 

Xa 

ta 

Qrs 

Vrs 

Urs 

/r 


Sirs 

^a,k 


yrrs 

T 

grs 

g 

K 

i^rs 

jrs 


node (index) set 

arc (index) set 

set of origin nodes; TZ G AT 

set of destination nodes; S C Af 

set of rontes connecting 0-D pair r — s; r G 71, s E S 

decision variable, independent of time 

decision step 

flow on arc a; x = (• • •, • • •) 

travel time on arc a;t = (■ • •, ta, ■ ■ ■) 

traffic demand between origin-pair rs; qrs = /o* 

average travel time for 0-D pair r — s 

travel time fnnction for variable demand for 0-D pair r — s 

flow on ronte k connecting 0-D pair r — s] = (/g^, • ■ •, fj^^, 

travel time on ronte k connecting 0-D pair r — s] = {■■■, c 


^■•■);c=(■•■,c-,••■) 


indicator variable: = 


1 if link a is on ronte k between 0-D pair r — s 
0 otherwise; 


the set of satisfying ([T]) 
the set of f; JF = 

the set of initially non-empty satisfying ([1]) and P^(r) > 0 for r = 0 
the set of f; ^ = Hrs 

the nnmber of initially non-empty rontes connecting 0-D pair r — s 
FIFO violation for flow on ronte k connecting 0-D pair r — s; 

J rs _ / jrs 

~ V ' ' ^k ’ 


|J ||2 2-norm of J(f), dehned in ffTSjl 


Table 1: Network notation system for static TAP 


Initial gness of f, e.g. iP-shortest routes for each 0-D pair 
for n = 1, 2, 3 ■ ■ ■ 

Find an equilibrium with Euler’s method 

If no shorter unused routes, the equilibrium is UE; 

Otherwise, perturb the equilibrium in the directions of shorter unused routes 
Table 2: A perturbation-based method for Ending stable UE 
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cumulative flow 



Figure 1: An example of 0-D First-In-First-Out violation 


Link 1 



tj=10(l + 0.15(|-r) 
t2 = 20(l + 0.15(^)^) 
t3 = 25(l + 0.15(^r) 

Xj + X2 + V3 = 10 


Figure 2: A static road network 
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Figure 3: Convergence of Euler’s method 



Figure 4: Solution trajectory of Euler’s method with r = 0.02 
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equilibria 

link flow 

link cost 

1 

10 

0 

0 

947.5000 

20 

25 

2 

0 

10 

0 

10 

137.1875 

25 

3 

0 

0 

10 

10 

20 

487.9630 

4 

4.0346 

5.9654 

0 

34.8405 

34.8405 

25 

5 

4.7864 

0 

5.2136 

59.2053 

20 

59.2053 

6 

0 

6.0762 

3.9238 

10 

35.9740 

35.9740 

7 

3.5833 

4.6451 

1.7716 

25.4560 

25.4560 

25.4560 


Table 3: All equilibria for network in Figure [2] 



Figure 5: Stability of UE and instability of PUE 














7^ 

5 

ICrs 

T 

Ar 

t 

[0,To] 

[0,T] 

N 

M 

At 

At/M 

tn 

u 

frir, t) 

9k{r,t) 

fris, t) 

9k{s,t) 

Prs{r,t) 

qrs{r,t) 

Prs{s,t) 

c"k{r,t) 

Vrs{r,t) 

Ht) 

K 

j.^rs 

IIJII2 


set of origin nodes; IZ d N 

set of destination nodes; S <Z Af 

set of routes connecting 0-D pair r — s; r E TZ, s E S 

independent decision variable 

decision step 

independent time variable 
assignment interval 
traffic simulation interval 
the number of assignment time steps 

the number of traffic simulation steps during each assignment time step 
assignment time step, At = Tq/N 
traffic simulation time step 

assignment time instants, tn = nAt for n = 0, ■ ■ ■ iV 

traffic simulation time instants, ti = iAt/M for i = 0, - ■ ■, MNT/Tq 

arrival flow at origin r at t taking route k connecting 0-D pair r — s; 

r*(r, f) = (• • ■, f), ■ ■ ■); f(r, f) = (• ■ ■, r"(r, f), • • ■) 

arrival flow-rate at origin r at t taking route k connecting 0-D pair r — s; 

g"(D t) = {■■■, gl^{r, t), ■ ■ ■); g(D ^) = (''', g''"(D 

arrival flow at destination s at t taking route k connecting 0-D pair r — s] 

arrival flow-rate at destination s at f taking route k connecting 0-D pair r — s] 

arrival flow at origin r between origin-pair r — s; 

arrival flow-rate at origin r between origin-pair r — s] 

arrival flow at destination s between origin-pair r — s; 

arrival flow-rate at destination s between origin-pair r — s; 

travel time on route k connecting 0-D pair r — s for a vehicle arriving at origin r at t; 
travel time between 0-D pair r — s for a vehicle arriving at origin r at t; 
the set of satisfying flTBl) 
the set of f(f); J^{t) = 

the number of initially non-empty routes connection 0-D pair r — s 

FIFO violation of flow from destination r at f on route k connecting 0-D pair r — s; 

r^{r,t) = J(r,f) = (• • •, v^"(r, f), • • ■); J = (■ • •, J(r, f„), • • ■) 

2-norm of J(f), dehned in (1^ 

Table 4: Network notation system for dynamic TAP 
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Initialize Tq, T, r, N, M, At, At 




Initial guesses of tn) V k, r, s, n 




for Ti = At, 2 At, 3 At ■ ■ ■ 




Use traffic ffow model to obtain 

ti) V k,r,s 

and i 

= l,---,MNT/To 

Compute travel time d'if{r,tn) and Vrs 

(U tn) 



Compute FIFO violation functions Jl‘^{r,tn) V k. 

r, s 


Use Euler’s method to solve gl,^{r,tn) 

and hud fj)' 

'(U tn) 

{W n = l,---,N) 


Table 5: An algorithm for finding a dynamic equilibrium 
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Link 2: Length=2 L 




Figure 7: A dynamic road network 



Figure 8: A dynamic UE for the network in Figure [7] 
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Figure 9: UE arrival flow-rates at the origin and travel times for the network in Figure [7] 



Figure 10: Convergence of solutions to a dynamic UE with different initial conditions with 
At = 0.05 for the network in Figured 


32 















